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Interior perfect fluid solutions for the Reissner-Nordstrom metric are studied on the basis of a new 
classification scheme. It specifies which two of the characteristics of the fluid are given functions and 
accordingly picks up one of the three main field equations, the other two being universal. General 
formulae are found for charged de Sitter solutions, the case of constant energy component of the 
energy-momentum tensor, the case of known pressure (including charged dust) and the case of 
linear equation of state. Explicit new global solutions, mainly in elementary functions, are given as 
illustrations. Known solutions are briefly reviewed and corrected. 
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The unique exterior metric for a spherically symmetric charged distribution of matter is the Reissner-Nordstrom 
^0 \ solution. Interior regular charged perfect fluid solutions are far from unique and have been studied by different authors. 
■ The case of vanishing pressure (charged dust (CD)) has received considerable attention. The general solution, in which 
the fluid density equals the norm of the invariant charge density, was presented in curvature coordinates by Bonnor 
0. The proof that this relation characterizes regular CD in equilibrium, i.e. in the general static case, was given 
later ^,|). In the spherically symmetric case another proof was proposed in Ref. Concrete CD solutions were 
studied in these coordinates PH. The generalization of the incompressible Schwarzschild sphere to the charged 
case with constant T® was also undertaken in a CD environment ||. Charged dust, however, has been investigated 
more frequently in isotropic coordinates, since these encompass the entire static case and allow to search for interior 
solutions to the more general Majumdar-Papapetrou electrovacuum fields [^,^|. In both coordinate systems there is a 
,_H ' simple functional relation between goo and the electrostatic potential. In isotropic coordinates there is one non-linear 
main equation which has been given several spherical pO|-|T3 | and spheroidal [ pl|^2|Jl^ ] solutions. One of them 
' coincides with the general static conformally flat CD solution . These CD clouds may be realized in practice by 
q-i a slight ionization of neutral hydrogen, although the necessary equilibrium is rather delicate. They have a number of 
interesting properties: their mass and radius may be arbitrary, very large redshifts are attainable, their exteriors can 
be made arbitrarily near to the exterior of an extreme charged black hole. In the spheroidal case the average density 
I* ■ can be arbitrarily large, while for any given mass the surface area can be arbitrarily small. When the junction radius 
ro shrinks to zero, many of their characteristics remain finite and non-trivial. One can even entertain the idea for a 
point-like classical model of electron were it not for the unrealistic ratio of charge e to mass m Q| . 

Recently, new static CD solutions were found, in particular with density which is constant or is concentrated on 
thin shells 16 . In the spherically symmetric case a relation has been established with solutions of the Sine-Gordon 
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The necessary condition for a quadratic Weyl-type relation has been derived also for perfect fluids with non- vanishing 
pressure However, in this case many other dependencies between the electrostatic and the gravitational 

potential are possible, even when combined with constant Tq |pl|| . 

The original Schwarzschild idea of constant density has been also tested in the charged case for a perfect fluid Q 
or for imperfect fluid with two different pressures [[l]. An electromagnetic mass model with vanishing density has 
been proposed in Ref. Unfortunately, the fluid has negative pressure (tension). Although the junction conditions 
do not require the vanishing of the density at the boundary, this is true for gaseous spheres. A model with such type 
of density was proposed both in the uncharged and the charged case [^5|,^0j . 

Another idea about the electromagnetic origin of the electron mass maintains that, due to vacuum polarization, 
its interior has the equation of state p + p = 0, where p is the density and p is the pressure. This leads to tension, 
easier junction conditions and realistic e and m p7|-p9[. It can be combined with a Weyl-type character of the field 
p0| . The experimental evidence that the electron's diameter is not larger than 10~ 16 cm, however, requires that the 
classical models should contain regions of negative density |5l],[52| ■ Probably an interior solution of the Kerr-Newman 
metric is more adequate in this respect. 

The presence of five unknown functions and just three essential field equations allows one to specify the metric and 
solve for the fluid characteristics [p3[. This is impossible in the uncharged case. Another approach is to electrify some 
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of the numerous uncharged solutions. This has been done for one of the Kuchowicz solutions |34] in Ref. |35(]. Two 
other papers [j36|]37f l build upon the Wyman-Adler solution (3^j3^]. Thus a charged solution is obtained, which has 
approximately linear equation of state when m/ro is small. In Rcf. |40) a generalization of the Klein- Tolman (KT) 
solution |^l|,^2| was performed, but the resulting fluid does not possess a linear equation of state. Recently, static 
uncharged stars with spatial geometry depending on a parameter (4^^^] have been generalized to the charged case 

The purpose of the present paper is to present a new and simple classification scheme for charged static spherically 
symmetric perfect fluid solutions. The calculations in each case are pushed as further as possible and general formulae 
are given in many instances. The known solutions are reviewed and compartmentalized according to this scheme in 
order to illustrate general ideas, without being exhaustive. New solutions are added where appropriate. The intention 
has been to stick to the simplest cases and remain in the realm of elementary if not algebraic functions. The junction 
conditions and other requirements for physically realistic models are discussed. The emphasis is, however, on the 
general picture, which appears unexpectedly rich and simpler than in the uncharged case. 

The metric of a static spherically symmetric spacetime in curvature coordinates reads 

ds 2 = e v dt 2 - e x dr 2 -r 2 dtf, (1) 

where dfl 2 is the metric on the two-sphere and v, A depend on r. The fluid and its gravitation are described by 
five functions depending on the radius: A, v, p,p and the charge function q which measures the charge within radius 
r. There are only three essential field equations, hence, two of the above characteristics must be given. We shall 
classify solutions according to this feature. For example, (y 1 A) is the case of given metric and the other three fluid 
characteristics are found from the equations. This does not mean that solutions are distributed among groups which 
do not overlap. Thus, (p, q) is a completely general case - any solution, after p and q are known, may be put into this 
class. The essence is that p and q 2 are given functions and there is control over them; they can be chosen regular, 
positive and comparatively simple. Then the other three functions are usually more complex and are not always 
physically realistic. 

In Sec. II the Einstein-Maxwell equations are organized into three main and two auxiliary ones. Two of the main 
equations are universal for all cases, the third one varies from case to case. Cases with given v have linear first-order 
differential equations for e _ V Cases with given A have linear second-order equations for e v l 2 or non-linear first order 
equations for p + p. These results hold also for fluids with a linear equation of state. The difficulties in the cases {p,p) 
and (p, q) are discussed, which prevent their analytical treatment. The junction conditions are given in general form 
and a reasonable set of physical requirements is included. 

In Sec. Ill the case (A, v) is briefly discussed. More attention is paid to the condition Tq = const which is a case of 
given A. The simplest algorithm for generating solutions is presented. It is based on a master function and involves 
simple operations. The same function is relevant also when v is given. 

In Sec. IV the case p + p = is solved completely. It represents a charged generalization of de Sitter spacetime. One 
new solution is worked out as an example. This case is generalized to include solutions with positive pressure and an 
explicit example is given. 

Sec.V deals with the case of given pressure which includes as a subcase charged dust solutions (they have vanishing 
pressure). Efforts are concentrated on the case (v, p) which is solved in general, A given as an explicit functional of 
the pressure. Conditions for regular density are elucidated. The well-known case of regular CD is used as an example 
to illustrate how the classification scheme works. Known solutions are reviewed both in curvature and isotropic 
coordinates. A new solution is given to demonstrate the features of a simple ansatz for v, used with variations in 
the next sections, too. The most general CD solution is derived. It is not regular at the centre but can form a 
'halo' around a typical solution with given pressure. An explicit example of such triple-layered model with realistic 
properties is applied. 

The case (v, q) is discussed in Sec. VI. It is the most direct generalization of the uncharged case. An integral 
representation for e~ A is given. A series of ansatze for v is introduced in general, the first two cases are worked out 
explicitly. A new solution is presented and some known results are corrected. 

In Sec. VII the important case (y, n) of a fluid with linear equation of state, involving the parameters n and po, and 
given v is investigated. Unlike the situation for uncharged fluid, the third main equation is linear and of first order. 
It allows a generic integral representation for e~ A which comprises also the cases (v,q) and (y,p). The structure of 
the density and its regularity are discussed which imply certain restrictions on v. The solutions are in some way 
electrified versions of the Einstein static universe (ESU) . A basic simplification technique leads to an infinite series of 
models with realistic physical properties. Two of them are studied in detail, including a model related to incoherent 
radiation. This kind of fluid appears also as a special case of the integral representation. Another special case is also 
soluble in elementary functions but has negative pressure and is closer to the cases studied in Sec. IV. Finally, a third 
special case, which is a twisted generalization of an obscure integrable uncharged case, is shown to be ill-defined in 
the presence of charge. 
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Sec. VIII is dedicated to the case {v,p), closely related to (l>, n). A solution in terms of hypergeometric functions 
is found when the density is constant. An algebraic ansatz for the density, depending on two parameters, leads to 
realistic solution except for q 2 which turns out to be negative. The process of fine-tuning, however, may give truly 
realistic solutions in the future. 

In Sec. IX the three cases (p, </), (A, p) and (A, q), which are closely related to each other, are discussed. It represents 
a critical review of Refs. Hj2^,|2^,^,g6| . A unified scheme is followed, which fully explores the connections between 
the different subcases. Several new solutions arc obtained which demonstrate how hard is to escape the curse of 
negative pressure, which plagues purely electromagnetic mass models. The solution of Wilson |23| is formulated as a 
hypergeometric function. 

In Sec.X the theme about a linear equation of state is raised once again, this time in the frame of case (A, n). First, 
a natural generalization of the KT solution is obtained for constant A, which is quite different from the solution in 
Ref. ^0|. Then the ansatz for A, which ensures that Tq is constant, is studied in detail. The general solution of the 
second order linear equation is expressed once more through the hypergeometric function and degenerate cases in 
elementary functions are sought for. Three different classes are studied. One reduces to ESU. The second has nice 
orthogonal polynomial solutions for e"/ 2 but with negative pressure, much like de Sitter solution which is contained 
as a subcase. This unpleasant feature is also shared by a trigonometric solution for the analog of incoherent radiation. 
The study of the non-linear first order equation produces even more bizarre results. 

Sec. XI contains some discussion and conclusions. In all sections old solutions, known to the author, are reviewed 
and classified. New solutions have been checked by computer for realistic properties. 



II. MAIN EQUATIONS AND CLASSIFICATION 

The Einstein-Maxwell equations are written as 



tf^«P+5 = Te-* + i(l-e->), (2) 



*P-5 = ^- A -3(l-e- A ). (3) 



q 2 _ x (v" AW V n i/-X' 



(4) 



where the prime means a derivative with respect to r and k — 87rG/c 4 . We shall use units where G = c = 1. The 
charge function is obtained by integrating the charge density a. We shall use, however, q as a primary object and 
then 

nc = %e^ 2 . (5) 
Spherical symmetry allows only a radial electric field with potential <f> given by 

F 01 =(j>' = -^e^. (6) 
Eq.(2) may be integrated by the introduction of the mass function 

M (r) = £ jf r 2 T°dr = £ jf [p + r 2 dr, (7) 

and gives 

, se -* = l-^. (8) 
r 



This can be rewritten as 



K P+^4 = — = ^C 1 - 2 -^)' ( 9 ) 



3 



and constitutes the first of our main equations. The second is obtained as a sum of Eqs. (2) and (3): 

K(p + p) = — (v> + \') = tv> --. (10) 
r r r 

The third main equation will not be Eq.(4) but another combination of Eqs.(2)-(4) which varies from case to case. 
One can transform Eqs. (2)-(4) into expressions for p, q, p 

, (v" AV v' 2 V 3i/ 1\ 1 
2g 2 _ _ A (v" X'v' v ,% X' v' 1 \ 1 



, , i/" AV i/' 2 5A' i/' 1 \ 1 

2«p = e- A - — + — - + — + - t +-o- 13 

i 2 4 4 4r 2r r 2 ' "- 2 



These equations may be written as linear first-order equations for z, suitable for the cases (v,q) ,(v,p) and (y,p). 
Introducing y = e v l 2 we have 

4q 2 

(rV + ry) z' = -2 (r 2 y" - ry' ~ y) z ~ 2y + -^-y, (14) 

(r 2 y' + hry) z' = -2 (r 2 y" - ry' + y) z + 2y - Anpr 2 y, (15) 

(r 2 ?/ + ry) z' = -2 (r 2 y" + 3rt/' + y) z + 2y + 4npr 2 y. (16) 

In the uncharged case the prescription of an equation of state makes the system of field equtions extremely difficult 
to solve. This is true even for the simplest realistic linear equation of state 

p = np-p , (17) 

where n is a parameter taking values in the interval [0, 1] for physically realistic solutions, while po is a positive 
constant, allowing the existence of a boundary of the fluid where p — 0. When po = we obtain the popular 7-law 
(with notation n = 7 — 1). In this case Eqs. (2)-(4) with q = lead to an Abel differential equation of the second kind 
independently from the approach or the coordinate system It is soluble in few simple cases. Almost all of the 
known few solutions of the more general Eq. (17) may be obtained by imposing a simple ansatz on A which makes 
the system overdetermined f47| . Therefore, it is surprising that in the more complex charged case fluids, satisfying 
Eq. (17), are subjected to a linear equation, similar to Eqs. (14)-(16). Plugging Eq. (17) into Eq. (10) yields 



npy 



Replacing this result in Eq. (15) leads to 



2 np a )y 

r V r 



(18) 



r^' + — ry ) z' = -2 r 2 y" + ——ry' + y)z + 2y f-r 2 y. (19) 

n + 1 / \ n + 1 / n + 1 

We call this case (v, n). It is a hybrid between (is, p) and (v,p) because fixing the equation of state, we fix the pressure 
in terms of the density. Eq. (17) is a privileged one, due to its linearity. Another realistic equation of state, namely 
the polytropic one, reads p = np 1+1 / k . Replaced in Eq. (10) it gives the relation 



, , 5n + 1 \ , _ ( 2 ", 3 - n / , \ , o 4k P° 2 



k( P + np 1+1 ' k ) =^-i (20) 
V J ry r 



This causes the appearance of radicals in Eq. (15) at best, and to non-integrable equations. 
Eqs. (14)-(16) and (19) may be written in the general form 
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gz' = f lZ + f , (21) 

whose quadrature is 

z = e F (C + H), (22) 

F= [ J±dr, H = [ e- F —dr. (23) 
J 9 J 9 

Here and in the following C will denote a generic integration constant. They may be written also as linear second-order 
differential equations for y, useful in the cases (A, q) , (A, p) , (A,p) and (A, n): 

2r 2 zy" + (r 2 z' ~ 2rz) y' + (rz' -2z + 2- y = 0, (24) 

2r 2 zy" + {fz 1 - 2rz) y' + {hrz 1 + 2z - 2 + Anpr 2 ) y = 0, (25) 

2r 2 zy" + (r 2 z' + 6rz) y' + (rz' + 2z - 2 - Anpr 2 ) y = 0, (26) 



2r 2 zy" + r 2 z' + 2 rz ) y' + rz' + 2z - 2 H y = 0. (27) 

V n + 1 / \n+l n + 1 J tf v ' 

The coefficient before the second derivative is one and the same in all cases. Eq. (24) is the generalization to the 
charged case [^|[|7| of the Wyman equation Obviously, the case n = — 1 is not covered by Eqs. (19), (27). 

One can find first-order differential equations also for the cases (A,*), based on the well-known Tohnan- 
Oppenheimer-Volkoff (TOV) equation generalized to the charged case Ssl. Let us first derive briefly TOV. 

Eq. (3) may be written as an expression of v' in terms of p, q and M 

, = 2M + upr 3 - q 2 /r 
V r{r-2M) ' 1 ' 

Now, let us take the following combination of equations (3)' + 2 (3) — - (4). After some cancellations it gives 

k P ' = - K (p + p)^ + ±2J-. (29) 

Inserting Eq. (28) into Eq. (29) yields the generalized TOV equation. We can trade q in Eqs. (28)- (29) for p and A 
by using Eq. (9). The result is a Riccati equation for Y = k (p + p) 

Y> = --e x Y 2 + -Y - + - & ir 2 M')' . (30) 
2 2 r r 4 

It marks another way to solve the cases (A, p) , (A, n = — 1). Its solution yields for the pressure up = Y — up. Then 
v is found from Eq. (10), while q is found from Eq. (9). Eq. (30) is a non-linear, but first-order companion of Eq. 
(15). Using Eq. (9) we find again a Riccati equation for the case (A, q) 

r x 9 A' 4q 2 8M' 2 , , „, 
From the definition of Y equations for the cases (X,p) and (A,n) follow 



Y> = - r -e'Y 2 + ( ± - t ) Y + 4«? + 4 (r 2 M>)' , (32) 
2 \ 2 r I r r 4 v ' 



Y' = - r -e>Y 2 + (*- , \ ) Y , iKP \ + A (r 2 M>)' . (33) 
2 \2 (n + l)rj (n + l)r r 4 V ; V ; 
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In all cases Riccati equations are obtained with the same coefficients before Y' and Y 2 . They may be written as 



Y' = --e x Y 2 + f 2 Y + f 3 , (34) 

where f% and f 3 are functions of r, differing from case to case. An arbitrary Riccati equation may be transformed 
into a linear second-order equation. The linearization of Eq. (34) is done by the change of variables u = e x ^ 2 y and 
reads 

- + A' + / 2 ) v! - T -e x hu = 0. (35) 
r J 2 

At this point f 3 contains A" which does not cancel. When a passage from u to y is done, A" cancels and we obtain 
exactly Eqs. (24)-(27). In this process we have exchanged Y, which is the sum of the pressure and the density, for y 
which is part of the metric. 

So far we have reformulated the original system of Eqs. (2)-(4) into Eqs. (9), (10) and a third equation, presented 
in many different forms, adapted to the cases of the proposed classification. We have briefly discussed the cases 
(A, v) , (A, *) , (v, *). The three remaining cases {p,p) , (p, q) and (p, q) are the most natural ones since one prescribes 
two of the fluid characteristics, hoping that the third one and the metric will be regular and reasonable. The case 
(p, q) is easily reduced to (A, q) because of Eq. (9). In the case (p,p), Y is also known and Eq. (30) becomes an 
equation for M 

2 , 9 ,\ i rW - M r 2 Y 2 , 4ko 

? ( ) + v^~m Y - w^m - Y -^-°- (36) 

This is an intricate, non-linear, second-order equation, which is not simpler than the TOV equation. It seems that it 
can be dealt with only numerically. 

The third main equation for the case (p, q) is obtained in the following way. Let us replace the density in the TOV 
equation (29) with its expression from Eq. (9). After some tedious manipulations, the following equation for M is 
found 

(M + go) M' = UM + f 5 , (37) 

K q Q 2 , 

°° = 2 pr -b (38) 

h = ^p'-^- K -r 2 p + ^ 21 (39) 

2/ 5 = -^(^-^)- 2 (^-^)(^4). (40) 

There is a standard procedure for the solution of such equations (4^^] ■ It consists of two changes of variables which 
bring them to the Abel equation of the second kind 

uj^ - u> = f (C) , (41) 

where lj = M + go while 

C = go + I Udr, (42) 



/(0 = ^xr. (43) 

9o + U 



Its integrable cases are few, depend on the shape of / (() and are tabulated in Rcf. 
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As a whole, the most attractive are the mixed cases (y, *) where one fluid characteristic and one metric component 
are specified. They lead to the most simple Eqs. (14)-(16) and (19). 

The five functions which describe the fluid together with its gravitational field should satisfy some physical require- 
ments. Eqs.(7)-(8) show that at the centre M (0) = and e A = 1. The density and pressure should be positive and 
monotonously decreasing towards the boundary. It is obvious that q 2 should be positive, too. The boundary ro of 
the fluid sphere is determined by the relation p (ro) = where a junction to the Reissner-Nordstrom (RN) metric 

x 2m e 2 .... 

e = er = 1 - — + ^2 ( 44 ) 

should be performed. The metric and v 1 must be continuous there. This leads to the expressions 



m = _ / rofA = Af (r ) q 2 (r ) 
ro V 2 J r 2r 2 



l_e"(l + ^l) =fim+«_^ ) (45) 



4 = l-e"(l + W) = ^^. (46) 
»o r 5 

The condition e = q (ro) follows from the vanishing of the pressure and vice versa. In fact, Eq. (10) gives at ro 

z = zv — nrop (ro) . (47) 

Replacing this in Eq. (9) and taking into account that e v (ro) = e~ A (ro), we get a proof of the assertion. 

Solutions, satisfying the above conditions are called in the following physically realistic. This is done in agreement 
with the majority of the reviewed papers, and their satisfaction is already rather non-trivial. Some other, more 
stringent requirements, are discussed in Refs. ||50|-|5^|. The most important of them is < ^ < 1, which means 
that the speed of sound is positive and causal, i.e., not greater than the speed of light. One has control over this 
characteristic in CD and models with linear equation of state. The case {y, q), which is a generalization of uncharged 
solutions, possesses their behavior for small q. In other cases the expressions for density and pressure may be so 
complicated that the fulfilment of this requirement is hard to estimate. We do not discuss these additional conditions 
in order to contain the present paper in a reasonable volume and prevent its dissociation into a set of separate 
publications. In many cases the examples have free parameters, or such can be added easily, so additional fine-tuning 
may be done. 

III. THE CASES (A, v) AND CONSTANT T ° 

In the case (A, v) Eqs. (11)-(13) should be used to find q,p and p. This is the simplest case but control over pressure 
and density is completely lost and one must proceed by trial and error. Krori and Barua ]33| have given the solution 

A = air 2 , (48) 

v = a 2 r 2 +a 3 , (49) 

where ai denote generic constants of a known solution. They are fixed by the junction conditions. The solution is 
non-singular and the positivity conditions are satisfied. 

Eqs.(7)-(8) show that the generalized Schwarzschild condition Tg = const determines A 

e~ A = l-ar 2 , (50) 

where a is a positive constant. We have the freedom to choose one more function. In Ref. || the condition p = was 
further imposed. We shall review this solution in Sec.V. In Ref. [2l|| some relations between cf> and v were utilized. 
They can be arbitrary in the perfect fluid case, so we do not base the classification on such criteria. One should use 
instead Eqs. (24)-(27) or Eqs. (30)-(33). The simplest expressions are obtained in the case (A, Y) when Eq. (30) is 
used to determine p 

r 2 Y 2 arY 

4np = -rY' - — + + 12a. (51) 

2(1 — ar z ) 1 — ar z 
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Then up — Y — up and 



g 2 = r 4 (3a - up) = — 

The function v is determined by a simple integration from Eq. (10) 

v 1 = — L-j (Y - 2a) . (53) 
1 — ar z 

This is an advantage over the second-order equations (24)-(27). Eqs. (50)-(53) solve the problem in a minimal 
algebraic way. Several positivity conditions follow for the master function Y and p: 3a > up, Y > 2a, Y > up , Y < 0. 

In order to illustrate how the above scheme works one may take any of the solutions elaborated in Ref. pl| , extract 
Y and check how the above equations and inequalities are satisfied. The simplest case has 

Y = 2a + ai(l-ar 2 ) 1/2 . (54) 

Of course, in any of the cases (A, *) we can study the subcase given by Eq. (50). This will be done in the following 
sections. 

A similar problem arises when v is given. Which function should be prescribed in addition in order to have the 
simplest algorithm for generating solutions? Again, Y is the best choice. This is seen by taking the difference of Eqs. 
(25), (26) or reformulating Eq. (10) 

yz — 2y' ' z + ryY = 0. (55) 
This equation is much simpler than any of Eqs. (14)-(16). 



Y' 



rY (Y - 2a) 
2(1 -ar 2 ) 



(52) 



IV. CHARGED DE SITTER SOLUTIONS AND THEIR GENERALIZATION 



In the uncharged case this is a special integrable case with linear equation of state (17) and n = — 1, po = which 
is equivalent to the de Sitter solution. The charged case is also completely solvable and falls in the case (A, Y = 0) 
according to our classification. The third main equation is Eq. (30) which becomes 

(r 2 M'\ 1 

^=H^ = 2^ (2M +rM >■ (56) 
Eq. (10) gives v — —A which is a feature also of the exterior RN solution. Eq. (9) yields 

, r 5 f M'\' r , . 
1 2 = - Y (jr) = 1 {2M'-tM"). (57) 

Thus, when M is given, all other unknowns follow from simple formulae. Eq. (57) allows also one to take q as a basis: 

2 pT 2 

K P = 2a -^~4 2_dr, (58) 
r Jo r 



M = -2 f f 2 r ^rdrdf + (59) 
Jo Jo rb 3 

Here ao is some positive constant and clearly q — r 2+£ q (r) where e > and g (0) = const. Eqs. (58)-(59) 
demonstrate the process of 'electrification' of de Sitter space. The bigger the charge function, the lower the density 
until some point ro is reached where p(ro) = and consequently the fluid sphere acquires a boundary. We already 
know that M > 0. Eqs. (56)-(57) show that p and q 2 are both positive only when M' > and 2M' > r \M"\. The 
equality holds at the boundary, where M" (ro) < 0. The charge and the mass of the solution follow from the junction 
conditions (45), (46) 

e 2 = 2r 2 M' (r ) , (60) 
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m = M (no) + r M' (r ) 



(61) 



Obviously, there is an abundance of solutions since M and its first two derivatives have to satisfy few simple 
inequalities. Three solutions are known in the literature. The Tq — const condition leads to q = in the interior 
and consequently to the de Sitter solution, which has constant density. One can introduce, however, a surface charge 
a ~ S (r — tq) which gives q — 6 (r — tq) e, p — po and 

z = l- |p r 2 . (62) 

This is the solution of Cohen and Cohen . Another solution j2q] has 

360 



M = ^—olr 3 (ba\ - 2r 2 ) , (63) 



q=^or 3 , (64) 
6 
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/c 2 



r A ) , (65) 



Z = 1 ~ + 90°° r ' (66) 

Here ctq enters the charge density a = a e~ x / 2 , while a\ is related to a Q = i^o\a\. This is one of the electromagnetic 
models of the electron. When Co — > we have p,p, q, A, v — ► due to ao — > 0. The limit is flat spacetime. We have 
shown, however, that ao is not obliged to vanish when q = so that, in general, the case (A, Y = 0) is an electric 
generalization of de Sitter spacetime. The natural generalization of flat spacetime is the CD solution, which will be 
discussed in the next section. 

A third solution has been found by Gautreau [^0|. It has 

2M = c 2 1 r-^-sm 2 — , (67) 
r a 2 

Z = 1 _ c l + ( C ^1^L]\ (68 ) 



c l „ s _a r_ 

CZ2 



(69) 



4> = ci H — ^— ^sin — , (70) 
r a2 

and q = —r 2 (j)'. Many other solutions are possible and we give one simple example, a variation of the solution in Ref. 

M = a(r 3 - r 4 ) , (71) 

up = 2a (3 - 5r) , (72) 

q 2 = 2ar 5 , (73) 

z = 1 - 2ar 2 + 2ar 3 . (74) 
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The junction at r = 3/5 gives m = arg, e 2 = 2ar^. M satisfies all necessary inequalities and the density is positive 
and decreasing. 

The case Y — is easy because Eq. (30) collapses into the simple relation (56). A generalization can be made by 
taking Eq.(56) as a basis. This constitutes the special case (A,/? (A)). Then q is given again by Eq. (57), while Eq. 
(30) becomes 

Y' = - r -e>Y* + ^Y. (75) 

This is a Bernoulli equation and, unlike the Riccati equation, it is readily soluble in quadratures. Its general solution 
gives an expression for the pressure 

KP= C+Ue^rdr - KP > (?6) 

where C _1 = Y (0). When C^oowe return to the previous case Y = 0, the trivial solution of Eq. (75). Fortunately, 
Eq. (10) can be integrated explicitly too and a closed expression is obtained for v 

e v ' 2 = A^e-^ 2 ^1 + -L jf e 3X/2 rdr^j , (77) 

1 r° 

i + 



1 r° 

± J e^rdr. (78) 



The second equation follows from the junction conditions. Eqs. (56),(57),(76)-(78) provide the generalization of the 
case Y = 0. Thus, every function M leads to two solutions: one with trivial Y and one with non-trivial Y, satisfying 
Eq. (75). The trivial solution has the disadvantage of negative pressure. Eq.(76) suggests that in the non-trivial case 
solutions with positive pressure may exist. 

Let us take Eq. (63) with somewhat different constants 

M = br 3 (1 - br 2 ) . (79) 

Then 

Kp = 36 (2 - 56x) , (80) 

q 2 = 56V, (81) 

up = — — + lbb 2 x - 66, (82) 

C(l- 26.x + 26 2 ; 2 ) 1/2 -i + f 

where x — r 2 and 6 > 0. When x < 2/56 the density is positive and decreasing. The pressure is positive at the centre 
when 1/46 < C < 5/126. Let us choose 6 = 0.01 and C = 40. Then p has a maximum at xq — 0.3 and a root at 
x = 2.5. This is a semi- realistic interior solution with e 2 = 56 2 Xq and mass given by Eq. (61). 

V. THE CASE OF GIVEN PRESSURE. CHARGED DUST 

In this section we discuss the cases (f,p) and (A,p). The pressure is considered a known positive function which 
decreases monotonously outwards and vanishes at the boundary of the fluid sphere. The simplest case is p — (the 
only reasonable case with constant pressure, unlike the cases with constant density). This represents charged dust. 
The case (u, p) is much easier since the third main equation (16) is linear and first-order. Something more, in Eq. 
(21) fi = —2g'. This allows to obtain a compact expression for z 



1 C 4k f r ( , rv'\ „ < 

9 H 9 H 9 / 1 + -tH e V 

(1 + ^) r 2 e-(l + ^) rV(l + f) io \ 2 J 



pdr. (83) 
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The knowledge of v allows to satisfy two of the junction conditions by choosing a function continuous at tq together 
with its derivative. Eqs.(9)-(10) provide an expression for q in terms of the known v,p and with z (v,p) from Eq. (83) 



o 2 

^ = 1 - z(l + rv') + K P r 2 . (84) 



The density follows from Eq. (10) 



Kp = — I v ) — up. (85) 



Its structure is clarified when z' is replaced by v 1 , using Eq. (83) 

2z 2 2zv' zv' zv" Anp , , 

Kp = —, tt ^ 1 t 7T H r ^—r ~ up. (86) 



;i+ r f i ) »• r(i+^) l+T 1 



2 



When r —* 0, z — > 1 and the third term on the right produces a pole unless i/ — > 2^or, which means z/" — > 2^o. 
Then rz^' — > and the poles in the first two terms cancel. The last two terms approach negative constants. In order 
to compensate them and have a positive density at the centre, vq must be a positive constant and the inequality 
8^o > 5np (0) should hold. Another consequence is that e v is an increasing function in the vicinity of the centre. 
Since e~ A is a decreasing function which starts from 1 and equals e v at tq, we have e v (0) < 1. Eq. (83) shows that z 
has a pole at r — unless C = 0. 

Now we can understand the physical meaning of the terms in z. The third term represents the contribution from 
the pressure to the metric. When p — 0, z is still non-trivial and represents the general CD solution. It has a pole at 
the beginning of the coordinates and should not be used there. When C = only the first term remains and this is the 
regular CD solution in curvature coordinates [Q. The intricate proofs that this is the most general regular solution 
are replaced here by the obvious condition C = 0. The second term in Eq. (83) may be induced in principle by 
the third if the lower limit of the integral is changed. The first term may be absorbed by the third if one uses the 
identity 



TV . , 

1 + — \e v r 



5(rV)', (87) 



2 J 2 

and makes the shift Kp — > Kp + gpr. The expression under the integral in Eq. (83) is bell-shaped, since it vanishes 
both at zero and tq and is positive in-between. 

Let us discuss in detail first the case of regular CD. It provides an excellent illustration of the classification scheme, 
adopted in the present paper. In the case (v,p = 0) Eqs. (5)-(6) and (83)-(86) give 



2(1 + ^)' 



(89) 



rv" + 2v' + \rv n 
73 ' 



Kp = K(7 = , (90) 



<t> = -e vl2 + ^o- (91) 
For the case (X,p = 0) we should use Eq. (88) to express everything via A 

v' = I (e x ' 2 - l) , (92) 

q = r(l- e- A / 2 ) , (93) 
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np = kg = A (V A/2 - e- A + ^A'e- A ) . (94) 
These formulae involve only the first derivative of A. The case (q,p — 0) is explicitly solvable, unlike the general case 

2 

? X = T^—2> (95) 
{r-q) 

M = q-£, (96) 

v 1 = -j^—y (97) 
r (r — q) 

2 

np = kg = — (r — q) q . (98) 

The condition q < r must hold. These formulae clearly demonstrate the electrification of flat-space, which is the 
trivial dust solution in the uncharged case. If we put p — in Eqs. (37)-(40) we still obtain a complicated Abel 
equation. One can check that M, given by Eq. (96), satisfies it. The most complicated is the (p,p = 0) case. If we 
introduce w via 



then Eqs. (92), (93) give 



e x = w- 2 , (99) 



= HlzHO , (100) 



q = r(l-w), (101) 

while Eq. (94) becomes an equation for w 

rww 1 = —id 2 + w — —r 2 p. (102) 

This is an Abel differential equation of the second kind like Eq. (36). The procedure for its solution was described 
briefly after Eq. (40) and brings it to the canonical form ]i^^9]] 

WW' -W = -^r 3 p, (103) 



e- x = 



where W = wr. The set of density profiles, leading to integrable W is very restricted. 

Several explicit dust solutions in curvature coordinates have been given. Efinger Q studied the simple case 
4/9 in the presence of a cosmological constant. When it is zero we have q — r/3 and singular e v — a\r, p = g = A/Qkt 2 , 
^10 = \\fa\r~ 1 l 2 . Florides illustrated his general discussion with a power-law solution q = a 2 r a3+3 , 



= {l-a 2 r a *+ 2 ) 2 , (104) 
(l - a 2 r a3+2 )^^ 1 , (105) 



p = o=a l> [L) V X ~. 00! 



where a 2 = 2(a*+3)b a3 ' F ma %; m ^ e f- lH tne case 01 constant Tq was discussed. It leads to Eq. (50) and Eqs. 
(92)- (94) give 
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v = 2 In 



0,4 



1 - (1 - a 2 r 2 ) 



1/2- 



(107) 



2 2 
q = r 



2-aV-2(l-aV) 1/2 



(108) 



K/9 = 4a 



2 r 



2 2\l/2 



1 - (1 - aV) 



(109) 



We use the known solutions just to illustrate the general ideas and formulae and do not elaborate on their physical 
significance or details of the junction conditions, which may be found in the original references. 
As mentioned in the introduction, most work on CD has been done in isotropic coordinates 



ds 2 = lJ- 2 dt 2 - U 2 (dr 2 



l dn 2 



where the simple relation 



±U 1 holds and the only equation to be solved is 
2, 



U" + -U' 



-U 6 p. 



(110) 



(111) 



When U is given, we have the mixed case (U, p = 0) and the density is readily determined from Eq. (111). The 
function U~ 2 = air 2 + 0,2 was used in Refs. |LG|L5|]. Another, more general function 



a 3 (r§ - r k ) 
kr k +1 



(112) 



was studied for k = 2 in Ref. [Q, for k — 4 in Ref. |jll| and for general k in Ref. |r|. Recently, the function 
U = a 4 MH W as studied @@. 

When p is given, it is convenient to transform Eq. (Ill) into 



u xx = — — — , 



up U 3 

Yx~ 4 



(113) 



where X = l/r. We may consider p{U) as a known function, chosen to simplify this equation. One possibility, 
leading to Bessel functions, is p = ae/U 2 (l^JT^]. Other choices lead to two-dimensional integrable models, like the 
sine-Gordon equation JTq| . The metric (110) describes also the regular static CD. Then the case of constant fluid 
density po is soluble in cn, one of the Jacobi elliptic functions, but the solution cannot be spherically symmetric. In 
the spherically symmetric case Eq. (113) becomes an Emden-Fowler equation, whose integrable cases are tabulated 
in Ref. p9| |. When the power of U is 3 like in Eq. (113), the integrable cases have X 3 or X 6 and not X 4 . We can 
make the transformation ip = —U/r, Z = U r — U/r and obtain 



77 ry K Po ,3 

ZZ4, - Z = — —ip . 



(114) 



This equation coincides exactly with the Abel equation (103) in curvature coordinates with constant density and is 
non-integrable too. 

Let us go back to curvature coordinates (1) which are more convenient when one wants to study also the case of 
non- vanishing pressure. From the considerations about the regularity of p, given after Eq. (86), it follows that the 
simplest function v would be v = vqt 2 . This choice coincides with Eq. (49) and leads to exponential functions and 
probably to the error function in Eq. (83). We shall choose a different option 



br 2 



(115) 



where < a < 1 and 6 > 0. This ansatz has been thoroughly studied in the uncharged case |30],[52|. The CD solution 
build on it has 



- br 
~2br2 



(116) 
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br 3 



2br 2 ' 



(117) 



up 



2b (a + br 2 ) (3a + 2br 2 ) 



(a + 2br 2 f 

Density is monotonously decreasing. The junction conditions (45)-(46) give 



brl 



a + 2br 2 



2 ' 



(m 



(119) 



m = r (l 

The condition z — e v at r supplies the relation 



2brl) . 



1 - 4a + VI + 8a 



' o 



(120) 



(121) 



which ensures the obligatory e = m. The r.h.s. is positive when a < 1, which is the case. 

The main disadvantage of regular CD solutions is that they possess a fixed ratio e/ra which is unrealistic, especially 
for classical electron models, and requires the extreme RN solution as an exterior. Eq. (83) tells that when the point 
r = is excluded, general CD solutions are possible. They have the following characteristics 



z=[l + 



1 + C- 



(122) 



up 



rv 

1 + - 



1 + C- 



z 
r 



(123) 



4_,_ >±=L( 1+C «; 

r 2 (l + ^Y V r2 



(124) 



h , = c ,/2 [rV 2 -4C(lW)e-"] 
2r(r 2 +Ce-") 1/2 



(125) 



When C = they reduce to Eqs. (88)-(91). What is the physical significance of such solutions? Let us consider a 
core of perfect fluid, occupying a ball with radius ro, with given v and -p. Then z is determined from Eq. (83) with 
C = 0: 



1 



i 



C(rY 



(126) 



C (r) = 4k (l+ r -^j e u r 3 pdr. (127) 

At the junction z is equivalent to Eq. (122) where C = Co = C (ro). An observer cannot understand whether the 
interior solution consists of perfect fluid or general CD since the only imprint left by the pressure distribution inside 
is a constant. In principle, the RN metric should be taken as an exterior, but there are no obstacles to take a general 
CD metric and postpone the junction to RN till another point n > r . Thus we obtain a tripledayered model with 
a perfect fluid core up to r (zone I), a halo of general CD from r to n (zone II) and the RN solution for r > n 
(zone III). The metric component z is given correspondingly by Eq. (126), Eq. (122) with Co and Eq. (44). For 
this purpose we choose a continuous, together with its derivative, function v in the region < r < n. At the first 
junction z is also continuous, while the pressure drops to zero. The density is finite. In zone II the pressure remains 
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zero, while the density continues to decrease. Finally, at r% the density also drops to zero (perhaps with a jump) and 
a RN solution follows till infinity, making the composite solution asymptotically flat. As seen from Eqs. (123), (124), 
the constant Co shifts e/m from 1, like a perfect fluid solution, occupying zones I and II would do. 

This picture will be backed by an explicit example with the ansatz (115). Since the details should not depend on 
the form of p (r) but only on the value of C, we shall discuss first zone II and its junction with zone III. In the interval 
[rg, r{\ we have 



(a + br 2 ) [C + r 2 (a + br 2 )} 
r 2 (a + br 2 ) 2 



while in general 



z 

up = — 



8br 



C (r) + 2r (a + 2br 2 ) 
r^a + 2br 2 C (r) + ar 2 + br 4 



(128) 



(129) 



For CD C = and one can easily show that the density is positive for positive Co- It decreases towards zero when 
t — ► oo . The charge function is given by Eqs. (115), (124). Its square is positive when 



Co < 



'o 



3br 2 



The total mass and charge are given by Eqs. (45), (46) 

m = r\ (l — a — 2br\\ 



(130) 



(131) 



e 2 = r\ (1 - a- ibr\) . (132) 

The junction condition —A = v at r\ gives 

C = m 2 -e 2 . (133) 

This relation clearly shows the effect of pressure in zone I on the mass-charge ratio, not altered by the CD halo in 
zone II. The result in zone I was found in Ref. ||. 

The mass and squared charge are positive when 1 — a > 36r 2 . The constant Co is positive when 

Abr\ (a + br\) >a 2 -a + br 2 1 . (134) 

A sufficient condition is 

\ < br 2 < \ (1 - a) , (135) 

which requires also a < 0.24. Let us take a = 0.01 and br\ — 0.3. Then Eq. (130) is satisfied when (jr^J > 0.8 

and utilizing the equality we get r\ — 4.013\/Co, ro = 3.588\/Co, b — 0.075/\/Co. This proves the existence of CD 
solutions with |e| jm ^ 1. They sustain the value of C, obtained in zone I from a perfect fluid solution with positive 
pressure. 

It seems that there are no solutions of the type (y,p > 0) in the literature. It must be stressed that this case is 
completely general, unlike the cases of constant Tq or with p + p = 0, discussed in the previous sections. Every perfect 
fluid solution may be reformulated as a (v,p) case. Let us proceed with the ansatz (115) which transforms Eq. (127) 
into 

C (r) = 4k f (a + 2br 2 ) r 3 pdr. (136) 
Jo 

We base the discussion on C (r) as a fundamental object, hence, it is convenient to select a more involved p leading 
to a simple C (r). Let us take 

Po-Pir 2 
a + 2br l 
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and then 



C(r) 



Po 



(138) 



(a + br 2 ) 
{a + 2br 2 ) 2 



r 2 (po 



W 2 ) 



br 2 



(139) 



(a + 2br 2 ) 



b- bp 



3 aPl 



(140) 



where z follows from Eq. (126) while q is calculated from the general expression (84). The density is given by Eq. 
(129). Here po and pi are positive constants, connected by po = Pir 2 ,. The ansatz (115) yields = b/a and the 
condition p (0) > becomes 86/5a > up (0) = po/a. This is weaker than a necessary condition for positive q 2 , b > po- 
At the boundary the density is always positive since C (r ) — and Eq. (129) coincides with the expression for 
general CD. We have also po — 3Co/rg and using the numerical values of the CD example, one can express po and p\ 
in terms of Co, in addition to b, ro and r%. It can be shown that for C > 0.06, p and q 2 are positive and the solution 
is physically realistic. Thus we have constructed explicitly a triple solution of the type PF-CD-RN. This ends the 
discussion of case (v, p). 

The case (A,p) relies on Eq. (26). Even the simplest constant z leads to special functions or non-integrable equations 
for y when p is chosen physically realistic. 



VI. THE CASE (v, Q) 

The cases (v, q) and (A, q) comprise the most direct generalizations of the numerous uncharged solutions. Setting 
q = one obtains one of them with the chosen ansatz for v or A. Like the previous section, the case {v. q) is much 
simpler. The third main equation is Eq. (14). Now f\ is not simply — 2g' but it can be represented as a linear 
combination of g',g and y which take care for y" 7 y' and y respectively. This is a general result which holds for any 
of Eqs. (14)-(16), (19). In this particular case we have 

fi = -29' + - ~ 4y, (141) 
r 

H = 2 J e»(l + ^ e- 2h r~ 3 - lj dr, (143) 

f v'dr 

h = jTTW <144) 

and z is given by Eq. (22). We shall use an ansatz, slightly more general than (115) 

e v = (a + br 2 )\ (145) 



where k is an integer. In the uncharged case k = 1 leads to Tolman's solution IV J42|. The case k = 2 was discussed 
first by Wyman |38| as an extension of Tolman's solution VI and later was studied in detail by Adler [p9[. Solution 



with k = 3 was given by Heintzmann §53j (see also Ref. |52|). The general class was studied by Korkina |54[| , but her 



only explicit solution was the one of Heintzmann. Later Durgapal [|5J| studied in detail the cases k = 1 — 5. All of 
them satisfy the physical criteria used in this paper, but some have irregular behaviour of the speed of sound dp/ dp 
and p/p. 

Going to the charged case, Eq. (145) and Eqs. (45), (46) give 
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e 2 = q 2 (x ) , 



(146) 



r 1 + (2k + 1) txq \Xq 

where r = b/a, x = r 2 . Let us consider first the case k = 1. Then 

1 + ra; 



m_ 1 + J J Kfe + l)r^ 



l + 2rx 



(1 + + 2xQ) , (148) 



Q= / ^dx. (149) 



ii ^ 3 



The pressure and the density are given by 



1 + 2rx (1 + 2ra;) 2 1 + 2rx V x 2 J (1 + 2rx) 2 

r + C + 3tCx + 2(1 + 3tx) Q q 2 . . 

1 + 2tx x z 



When q = these expressions coincide with the pressure and the density from Refs. [}42|j55| l . We have Kp(Q) — 
3 (r — C). Let us choose q 2 = K 2 x 3 . Then the condition p (xq) — yields an expression for C 

c = _r+ (3 + 87*0)^*0 (152) 
1 + 3ra;o 

which means that C is negative and consequently p (0) is positive. The condition z = e v at xo defines a as 

_ 1 + Cx + 2AT 2 xg 



1 + 2tso 



(153) 



Together with Eq. (146) which becomes e 2 = K 2 Xq and Eq. (147) for k = 1 we can express c, a, e, m through xq, t, K. 
Let us consider next the case k = 2. Then Eqs. (142)-(145) give 

e F = X -^, (154) 

(a + 3te) 2/3 

z = l + Ce F + e F ( 2{a + bx \j dx . (155) 
7 (a + 3&x) 1/3 x 3 

When g = this is the result of Adler who used the simplifying assumption e F H = 1 instead of the ansatz (145). 
When q 2 = K 2 x 3 we obtain 

z=l+ - + ^^ 2 x+|i, 2 x 2 . (156) 

(a + 3bx) 2/3 56 5 

This is the corrected result of Nduka who wrote the field equations with an error in the sign of q 2 . In all his 
results K 2 should be replaced by —K 2 and the conclusions correspondingly altered. 

In another paper Singh and Yadav |57) simplified Eqs. (14), (21) by demanding fx = —fo- This leads to the 
following equation for y 

2a 2 

r 2 y" - ry' - -\y = 0. (157) 

Letting q = Kr they obtain the Euler equation. It has three types of solutions, all of which are singular at r = 
since y ~ r ai . If K = the non-singular solution of Adler is again recovered. 
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VII. THE CASE {y, N) 



This is the case with given v and the perfect fluid satisfying the linear equation of state Eq. (17). Realistic solutions 
have < n < 1 which means constant speed of sound and causal behavior. Surprisingly, the charged case is much 
easier than the uncharged one. The third main equation is Eq. (19) with n ^ — 1. An expression for /i, analogous to 
Eq. (141), may be written 



Then F and H read 



fi = -2g' + 



I6ng 



8n (9n + 1) 
—V- 



(n+l)r (n + 1) 



-(a+l) p /3hi 



- v (A 



rv' 
2 



(158) 



(159) 



H = 2 J e"(A+ r -^j r a e-^ (l - |jyr 2 ) dr, 



where 



hi 



~ J A+^L' 



(160) 



(161) 



a 



1 - 3n 
5n+ 1' 



(162) 



A = 



5n + l 
n+1' 



(163) 



4n (9n + l) 
(n + l)(5n+l)' 



(164) 



We shall study mainly the range < n < oo and occasionally some negative n. The case n = —1/5 demands special 
treatment to be done at the end of the section. The case n = — 1 was discussed in Sec. IV. Each of the coefficients 
a, (3 and A has the same limits when n — ► ±oo. When < n < oo their ranges are —3/5 < a < 1, 1<A<5, 
< < 36/5. Expressions (159), (160) are generic for the cases (v, q) , (u, n) , (u, p). The case (v, q), discussed in the 
previous section, is obtained by putting a = —3, = 2, A = 1, changing the sign of H and an obvious replacement 
for q instead of po- The case (v, p) will be studied in the next section. 

In the present case Eq. (162) shows that a + 1 is always positive for n € [0, oo]. Therefore, z will have a pole unless 
C = 0. Hence 



z = e F H. 

After z is found, Eqs. (9), (10), (17) allow to extract the fluid components from the metric 

1 



Kp - 



n+1 







np - 









(165) 



(166) 



up 







[t('-t) 


- npo 







(167) 



g 2 r 2 

n+1 



(168) 
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Since po > 0, a necessary condition for positive pressure is v' > z'/z and this leads to p > It is more convenient 
to express p and q from Eqs. (17), (84) respectively and lay all the difficulty in the calculations on the density. 
Plugging Eqs. (159), (160), (165) into Eq. (166) we obtain 



(n + 1) Kp = np + 



4np 



z {rv') 



z 
A 



2- 




(169) 



which is an analog of Eq. (86) and leads to similar conclusions for v. The poles in the two terms in the square brackets 
cancel when r — > because z — ► 1 and v' — » 2v r. This condition makes p a well-defined function. We shall use for v 
an ansatz which generalizes those in Eqs. (115), (145) 



e" = (a + &r")\ 
where s > 2 is not necessarily an integer. We obtain for z 



4 / (o + bxf- 1 [2aA + {2A + ks) bx} 1 (l - ^ff x 2 / s ) x^^dx 



sx^r- (a + bx) k 2 [2aA + (2A + ks) bx] 



1+7 



(170) 



(171) 



7 



2k(3 
2A + ks' 



(172) 



There are several ways to simplify the expression under the integral. One of them is to set 7 (which is not an integer, 
in general) equal to zero for any n choosing k and s appropriately. From Eq. (172) the following relation arises 



(36fc - 25) n 2 + (4k - 10) n - 1 



2 k(n+l)(5n+l) 
The R.H.S. should be > 1 which yields the condition n > 

k + 5 + A^k (2k + 1) 



= 



31fc - 25 



(173) 



(174) 



The first few values are n\ » 2.15, «2 ~ 0.53, 713 w 0.39, 714 = 1/3. The first is beyond the physical range, the fourth 
is an exact number. One cannot probe the region n < ps 0.21 with this method. 

Let us replace now z from Eq. (171) (with 7 = 0) into Eq. (169) and calculate p(0). When s > 2 the last two 
terms do not contribute, while from z only the term proportional to a k ^pQ counts. The result is very simple 



p(0) = 



3po 
3n + 1" 



(175) 
0U52I. 



The equation of state at the centre reads p (0) + 3p (0) = which is exactly the equation of state for the ESU [ 
The pressure is negative at the centre, which is unrealistic. 

In this connection it is worth to point out another relation to ESU even when 7 does not vanish. Let us take in v 
and z the limit 6 = 0. Then 



(176) 



KPq 2 

z = 1 -3^TI r ' (177) 

p is given by Eq. (175) and q — 0. Now we obtain ESU for any r independently of s,k,n and po. Hence, the results 
in this section may be considered as a specific 'electrification' of ESU like the generalization of de Sitter solution in 
Sec. IV and of flat spacetime in Sec.V. 

Going back to the case 7 = 0, we reach the conclusion that s = 2 is a necessary condition for realistic solutions. The 
overlapping bands of solutions with n > n% collapse to an infinite set of discrete values n^, k > 2. Then the ansatz 
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(170) coincides with (145) and we shall use the corresponding notation. A series of models is obtained, parameterized 
by k and n is given by Eq.(174). They have 



a k (1 + tf , 



(178) 



/(l + ^-^-^l-^jcft 
(l + t) k - 2 [A+{A + k)t] 



where p = kpq/t. Eq. (169) becomes 



ifx(l + t) 



up 

{U + 1} ~ =M+ ( n +l)[A+(A + k)t] ' t 



l + t 



2kz 
l + t 



2-/3 



l+t 



A + (A + k)t 



A A+(A + k)t_ 
Ikz 

(l + t) [A + (A + k) t] 



(179) 



(180) 



Now additional terms contribute to p (0) and p (0) and they can be made positive by a proper choice of r and po- The 
pressure is given by Eq. (17) and the charge function by Eq. (84) which becomes 



T<f = t 



1 + (2k + l)t up' 
l + t Z + ~ 



Let us discuss the particular case k = 2, ri2 = 0.53. Then 

_ 1.19 + 0.35* - 0.46^* - 0.27^< 2 
Z ~ 1.19 + 2.19* 

and it can be shown that p (0) > when p < 6.34. Let us take p — 1. Then 

k P _ -4.05 -6.1* +3+^ (4.57 + 6.57*) 
~ ~ 6.87+ 12.64* ' 



(181) 



(182) 



(183) 



rq 



til 



1 + 5* kp 
z + —t 

l+t T 



(184) 



The pressure is monotonously decreasing and vanishes at *o = 0.66. The charge function reaches a maximum before 
*o and is still positive at *o This indicates that a from positive turns to negative. As in Eq. (153) the condition 
z ( r o) = e 11 ( r o) determines a. The charge and the mass are determined by Eqs. (45), (46). All constants may be 
written as functions of ro: a — 0.37, b — 0.24/tq, m = 0.32ro, e = 0.49ro. The ratio e/m > 1. 

Another interesting case is k = 4 with = 1/3 exactly. When po — and q = this is the equation of state 
of pure incoherent radiation and the case is not integrablc fj^j . This specific charged version, which is not the only 
possible one, however, is integrable and assuming p = 2/3 we have 



1+1* 



I* 3 



A* 4 



(i + *r (1 + 3*) 



(185) 



up _ 5 + 9* 3 (3 + 7*) 2 | + 3* + f * 2 + i* 3 2 

V ~ 12 (1 + 3*) + 2 (1 + t) (1 + it) ~ 4 (1 + tf (1 + 3*) ~ 3 ' ( ' 

^..(x.^",^,). (18 r, 

The pressure is positive, monotonously decreasing and vanishes at to — 0.54. The charge density has a negative part, 
like the previous case. The junction conditions yield a = 0.3, b — 0.16/r 2 ,, m — 0.5ro, e = 0.45ro. This time e/m < 1. 
Probably, all models in the series have a range of p where they are physically realistic. 
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The method used to simplify Eq. (160), after all, resembles the Korkina-Durgapal method, which was extended to 
the charged case (v, q) in the previous section. With s = 2 the ansatz (170) is equivalent to (145) and, effectively, 
the non-integer 7 was exchanged for an integer k in Eq. (171), which allows integration in rational functions. There 
are other ways to simplify Eq. (171), even when k = 1. One of them is to take again n = 1/3. Then a = 0, A = 2, 
j3 = 3/2, 7=1/2. We have 

e" = a(l+t), (188) 



(l+t)/(2 + 3t) 1/2 (l-§ M t) t-^dt 



tV2 (2 + 3i) 



3/2 



The integral is expressed in elementary functions 
l + t 



z = 



t 1 / 2 (2 + 3t) 3/2 



(M - 4 + 3m*) V(2 + 3i)t + ^T^^ ln ( 3 * + 1 + V^t (2 + 3t)) 



(189) 



(190) 



where the integration constant was chosen so that z (0) = 1. The pressure and the charge are given by Eq. (17), (169) 
and (181) 



up _ 3(l + 2t) 9z 
rM + 



1 (z l + t 



4(2 + 3*)^" 4(2 + 3t) 2t \2 2 + 3t J ' 



(191) 



rq 2 = t 1 - 



3t 



l + t 



(192) 



The pressure is well-behaved for any /i and vanishes at to- However, q 2 also has a zero, but at ti, and becomes 
negative in a region where t > t\. Unfortunately, for most values of /i the inequality t\ < to holds, meaning that 
q 2 < in some part of the interior. This is grossly unrealistic. We have been unable to find by computer simulations 
any realistic solution. 

Another case with k = 1 which leads to simple functions can be guessed from Eq. (164). When n = —1/9, /? = 0. 
We have also a = 3, A = 1/2, 7 = 1. Although the pressure is negative, it is worth being discussed. Eqs. (188), (192) 
still hold. The integral in Eq. (171) is elementary and gives 



l + t 
(l + 3t) 2 



1 



„ 3 \ 27 2 



(193) 



Eq. (169) leads to p and p similar in form, since Eq. (17) is satisfied. Like usual, we give only the formula for the 
pressure 



up 

T 



l+t 



2(1 + zty 



3/27 
2+-M+ (3 + -H / 



9(1 + 2t) 3z 
4(l + 3t) M ~ 2(1 + 3*)' 



The values at the centre arc 



On the other hand 



f (0) = -i(l + 3 M )M, ^(0) = ?(1 + M ). 



np . . 1 9 
V (0O) = 6 - 8 M ' 



(194) 



(195) 



(196) 



Thus, when < fi < 4/27 the pressure is monotonously increasing, crosses the t axis and forms a boundary of fluid 
with tension inside. Let us take the 7-law with n = 0. The charge function q 2 is positive and increasing. The junction 
conditions give a — 0.55, b = O.II/tq, m = 0.23r , e = 0.35r where to = 0.2. This solution resembles the classical 
models of the electron in Sec. IV. 

Finally, some words should be said about the case n = —1/5, which is really very special. Eqs. (159), (160) are 
replaced by 
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H =\j eVr6 ( 2 " ^ Po r 2 ) e 2h2 dr, (198) 

f dr . , 

h2 = 2j-^. (199) 

Let us choose the simple ansatz (115) or (188) for v, written as 

e v = a + bx = a (1 + tx) . (200) 



Then Eqs. (197), (198) become 



b l x° 



where X = 1/x. The integrals can be expressed through the Ei (X) function. For example 

5 / \ — fe+1 fi— fe I \— 5 

«A (-r) x b (-t) 

(6 - 1) ... (6 - fe) 5! V rx 

Then z, which is given again by Eq. (165), will have poles for x = 0, not compensated by Ei. Thus, alreadyz is ill- 



X~ e e- x / T dX = e" 1 /™ y i iL- x ' + LzLi Ei ( - — ) . (203) 



defined, although the case is formally soluble. In addition n < 0. Whenpo = and q = 0, this case is soluble 46 1. It is 



and is best expressed in isotropic 



connected by a Buchdahl transformation to the de Sitter solution with n = — 1 
coordinates. This transformation, however, does not work for electrostatic fields and probably is not applicable to the 
present case of perfect fluid plus an electrostatic field. Nevertheless, some memories from the uncharged case seem to 
have been preserved in a rather twisted form in the presence of charge. 

The only paper which comes near to the results of the present section seems to be Ref. [g6| which belongs to the 
(y, q) case and is mentioned in Sec. VI (see Eq. (156)). When m/ro is small, this solution approximately satisfies Eq. 
(17) with n = 3/10, value rather close to 1/3. The connection to our solutions is not clear. 

A final remark is in order. It may be tempting to take the limit e — > in one of the solutions obtained here and 
try to find analytic expressions for the non-integrable at present cases of 7-law uncharged solutions p6| . However, 
for the ansatz (145) used in this section, this leads to ro — » and to flat spacetime. These solutions do not allow 
point-like models and this seems to be property only of the CD solutions, discussed in the introduction. In the case 
of general u, it seems that the above limit will produce charged solutions with zero total charge. This is suggested 
by the behavior of a in some of the models, namely, its change of sign in the interior. The proper non-charged limit 
seems to be q (r) = 0. Then, however, Eq. (9) relates p to A and we obtain the field equations of the uncharged 
fluid. Hence, nothing is gained in comparison with the three approaches described in Ref. [Q. They lead to Abel 
differential equations of the second kind with few integrable cases. 



VIII. THE CASE (y, p) 



This case arises as a special limit of (y, n) when n — > 00, po — > 00. One sees from Eq. (17) that in this limit 
P "~ *• Po = Po/ n which can be made a finite constant. More generally, the third main equation (15) is obtained from 
Eq. (19) plus the change ^_ _> p ( r ). Eqs. (162)-(164) give a = -3/5, A = 5, = 36/5. The constant C in Eq. 
(22) is zero again to prevent a pole at the centre. Eqs. (158)-(161), (165) yield 

h = ~2g' + ^ - 72y, (204) 
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o — 

2e 5 



5+ !^L) r -3/5 e -fft3( 1 _ 2K/9r 2) d? . 



(205) 



v'dr 



(206) 



The parameters n and po have been exchanged effectively for the density function and in this respect Eq. (205) 
resembles Eqs. (142), (143) from the [v, q) case. When p is constant we obtain a generalization of the Schwarzschild 
interior solution ^8|. It is much more complicated than the other generalization with constant Tq , which imposes 
just the simple ansatz (50) on the metric and was discussed in Sec. III. As seen from Eq. (205), the case p = is 
non-trivial and leads to purely electromagnetic mass models, because the mass function and T® receive contribution 
only from the charge. Formally, it is an analog of the 7-law equation of state (po = 0) from the previous section. 
Let us use again the ansatz (170) for v. Now Eq. (172) reads 



7 = 1 



12k 



50 + 5fcs' 



which is nullified when 



_ 72fc - 50 
S ~ 5k ' 

This expression always gives s > 2 for k > 1. The analog of Eq. (171) with 7 = is 

4 / (a + br s ) k ^ (l - 2npr 2 ) r -3 / 5 dr 
r 2 / 5 (a + br s ) k ~ 2 [10a + (10 + ks) br s ] ' 

Now p is given and p is found from Eq. (10) 



Kp - 



up. 



Let us calculate p(Q). When s > 2, just two terms contribute and the result is simple 



(207) 



(208) 



(209) 



(210) 



p(0) = --p(0). 



(211) 



This is similar to the relation following from Eq. (175). The pressure and the density cannot be both positive at the 
centre, hence, we must take s — 2. Then Eq. (208) shows that k is not an integer and we gain nothing, trying to 
simplify z. The method which lead to an infinite series of realistic solutions in the (y 1 n) case seems rather useless 
here. Therefore, let us put k = 1 in Eq. (171). Then 7 = — 1/5 and it becomes 



1 + t 



(1 - dt 
iV5 (5 + 6i) 4/5 J i 4 / 5 (5 + 6i) 1/5 ' 

When p = po the integrals bring forth the hypergeometric function 



l + t 



4/5 



F 



1 1 6 

5'5' 5 ; 



6f 
T 



ap °-tF 1 1 6 11 



3r 



5'5' 5 



Of 
T 



The charge is given by Eq. (192), while the pressure follows from Eq. (210) written as 

2z 



Kp 
T 



l+t 



2z t 



Kp 
T 



(212) 



(213) 



(214) 



Searching for an elementary solution, one may consider some sophisticated density profile, leading to a simple 
integral in Eq. (212). For example 



Kp 

T 



1 

2f 



1/5 



(1 -b 1 t + b 2 t 2 



(215) 
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leads to 



(i + ¥T 



6 *+lI* )' ( 216 ) 



where 61,62 are positive constants. We have p (0) > when 61 > 6/25 and p is decreasing when 62 > 4|p + 
Taking b± = 5, 62 = 2 a solution with realistic z,p and p is obtained, having a boundary at to — 0.52. Unfortunately, 
q 2 < in the interior, although for larger t it becomes positive. Perhaps, further fine tuning of the density may lead 
to physically acceptable solutions. 



IX. THE CASES (p,Q), (A,p) AND (A, Q) 

Having in mind the results in the previous section, it is no wonder that the case {p, q) , where there is control over 
the evasive charge function, is much more popular in the literature. Eq. (9) easily relates this case to (A, p) and (A, q), 
and Eq. (24) or (25) may be used as a third main equation. Eqs. (7), (8) define z in terms of p and q 

z = l-~ f pr 2 dr - - f trdr. (217) 
r Jo r J r z 

Eq. (9) shows that p is regular when z — 1 ~ r 2 for small r. We accept the ansatz 

z = 1 - ar 2 - br 4 . (218) 

It follows from Eq. (217) when q = Kr l with I — 2, 3 and 

p = Po - Pl r 2 , (219) 

where po,pi are non-negative constants.. When I = 2 or I = 3 we have respectively 

z = 1 - - (ko + K 2 )x + — x 2 , (220) 
3 5 



The total charge and mass are given by 



1^ x+ l -( Kpi -K 2 )x\ (221) 



Kx 3 /2 , (222) 



m = - [a+ (b + K 2 )x ] x 3 Q /2 . (223) 

When pi =/= this is a model of a gaseous sphere, both density and pressure vanish at the boundary |p5| , p6| . We have 
po = piXq. When pi =0, po this is a generalization of the incompressible fluid sphere for I = 3 or I = 2 |p3| . 
Finally, when P \ — po = we come to a model with electromagnetic mass |J. The case 6 7^ appears to be simpler, 
due to a peculiarity in Eq. (24) and leads to elementary functions. Eq. (50) demonstrates that the other case, 6 = 0, 
belongs to the family of solutions with constant Tq . Eq. (24) becomes in terms of x 

(1 - ax - bx 2 ) y xx - (bx + |) y x - Q + y = 0. (224) 

The coefficient before y is constant when I = 3 (or K — 0). Let us study this case first. Eq. (221) gives 




(225) 



A change of variables brings Eq. (224) to (see Eq.2. 1.164 from Ref. [|T|) 
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Vii + dy = 0, 



(226) 



Jo V 1 — ax — bx 2 



This integral has several different expressions according to the signs of b and a 2 + 46, but many properties of the 
solutions depend only on the facts that £ (0) = 0, £ and £ x are positive. Eq. (226) is easily solved and the last function 
to be determined, the pressure, is given by Eq. (10), suitably rewritten as 



up = Vi— - 2z x - up. (228) 



At the junction 



y{x ) = y/z(x ). (229) 
Then from Eq. (228) follows that the pressure vanishes when 

2% {x Q ) = z x (x ) (230) 
is satisfied. The charge is given by Eq. (222) while the mass in Eq. (223) becomes 



1 

m = — 
2 



K P0 , /VrX2 \ ^0 



_ , {QK 2 _ Kpi) ^ x * \ (231) 



When p\ = the mass is positive. When pi ^ this assertion still holds because then p±xo = pq. The positivity of 
the mass does not impose any conditions on the constants po, xq and K which determine the solution. 
Inserting z and p in Eq. (228) yields 

Kp = 4Vi— - \npo + - (3/cpi + 2K 2 ) x. (232) 
y 3 5 v ' 

Then 

«p(0) = 4^(0)-^«po (233) 

y 3 

should be positive. To be more concrete, let us introduce the three classes of solutions of Eq. (226) for d = 0, d < 
and d > respectively: 

V = Ci + C 2 t (234) 
y = C ie -^ + C 2 e^, (235) 

y = C 1 sin Vd£ + C 2 cos (236) 

The condition y > and Eq. (233) lead to C\ > 0, C 2 > when d > and to C 2 > 0, C 2 > C\ when d < 0. 
Therefore, e v increases with r until it meets z which decreases. We have y (0) < 1. 

1) Case d = 0. Eq. (225) means npi = UK 2 , while Eq. (221) gives b = -2K 2 . We also have 

and the solution is determined by two constants, p an d if- Eqs. (229), (230) result in 

C 2 = 4«Po, (238) 

DO 

Ci = ~ii K Po^ + V^, (239) 
66 
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Zq = 1 — —V , V 



Kpo 



3 ' n\K\ 

and y is given by Eq. (234). C\ is real when v < \/3/5. The charge is given by Eq. (222) and the mass is 

4 3/2 

m = — kpqXq 1 , 

leading to the ratio |e| jm = 3/4u > vT5/4. Due to the negative b Eq.(227) reads 



1 



:hl 



2bx 



In Ref. |26J| another expression was used, true only when 

a 2 +4b = —K 2 ( v 2 
9 



1 



72 
121 



:ln 



1 - 



< 0. 



This inequality, however, does not hold for v close to y/3/5. Plugging Eq. (242) into Eq. (239) results in 

lit; 



6V2 



In 



6\/2 



yfzo - In 



1 - 



6%/2 



The pressure reads 



up 



-Kpo + 7 A X. 



Ci + C 2 £ 3 



(240) 



(241) 



(242) 



(243) 



(244) 



(245) 



The condition p(0) > is fulfilled when C\ < 2/11. This holds when v varies near \/3/5 = 0.77. Thus, there are 
solutions with positive pressure. 

2) Case d < 0. This means repi < HX 2 and the limiting case pi — is possible. We shall discuss first the general 
case. Eqs. (221), (222), (231), (235) hold. Conditions (229), (230) fix the constants 



(246) 



° 2 - r 



Zx (go) 



The pressure is given by Eq. (232) 



up = A\/—dz 
Positivity at the centre demands 



77 ^ * ~ - Wo + - (3k P i + IK 2 ) x. 



C 2 - Cl Kpo 



C* 2 + Ci 12V^d 
A necessary condition is C 2 > Ci, satisfied when 

«Po 2 , 2 



^ (a?o) 



3 5 



(X 2 - Kpi) X 



is positive, which leads to npi > 6-ftT 2 . Introducing the variable 



A = 



K 2 ' 



(247) 



(248) 



(249) 



(250) 



(251) 



we obtain for its range 6 < A < 11. Then b < and £ is given again by Eq. (242). Eqs. (246), (247) become 
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Obviously, C2 > and C2 + C\ > 0, as required. Inequality (249) is still very complicated. It becomes simpler if 
d = b but this gives A = 3 which is out of range. Another way to simplify it is to seek solution with C\ — 0. Eq. 
(252) allows to rewrite this condition as 

_ 9 (11 -A) 

07 1 OA A2> 



27 + 9A- A 2 



while Eq. (249) becomes 



Plugging here the previous formula transforms Eq. (255) into 

(A + 2) (A -6) < 0, (256) 

satisfied when — 2 < A < 6. This interval is adjacent to the allowed interval for A, hence, this solution also possesses 
negative pressure. We have been unable to find a positive pressure solution. 

Let us discuss the subcase p\ — 0, i.e. p = po- Now Eq. (250) shows that z x (xq) < 0, therefore, C2 < C\ and 
either Eq. (249) is not satisfied or y (0) < 0, which is completely unrealistic. The pressure is negative in this subcase. 
Eq. (222) holds as it is, while Eq. (223) reads 

i K p + jU^o) x o 2 - (257) 

We have b = K 2 /5 > and consequently a 2 + 46 > 0. Then Eq. (227) reads 

1 a + 2bx 1 . a 

— = arcsin = arcsin — ■ 

Vb Va 2 + Ab Vb Va 2 + Ab 



1 . a + 2bx 1 . a 
£ = —= arcsin -= arcsin — (258) 



In Ref. |26j the expression 



1 a + 2bx , . 

i = — — arcsinh (259) 



f-b V-a 2 -4b 

which holds only when b < and a 2 + Ab < was used in all cases and this is incorrect. Formula (258) allows to make 
connection with the results of Ref. ||, where in addition po = 0, so that the density vanishes and the mass arises 
entirely from the electrostatic field energy. Eq. (250) still shows that the pressure is negative. Eqs. (222) and (257) 
give the ratio 

M = TOT-- ( 26 °) 

m 6\K\xa 

Now a = and 

z = 1 - \ r2x2 - ( 261 ) 

It is positive at the boundary when \K\ xq < y/E and this means |e|/m > v5/3. The second term in Eq. (258) 
vanishes and combining it with Eq. (235) yields 

y = C\ cxp I — arcsm Vbx J + C 2 exp arcsin Vbx J . (262) 
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The passage to arccos scales and interchanges their places, giving Eq. (4.6) from Ref. |(|. 

3) Case d > . Now we have np\ > UK 2 which permits the limit K = 0. The two cases do not differ essentially. 
Eq. (236) gives y. The constants read 

Ci = sin Vd£o H — - ^7=^ cos Vd£o, (263) 
2v d 

C* 2 - V^c° s Vd£o - sin Vd£ - (264) 

2y d 

The pressure is given by a formula, similar to Eq. (248). Positivity at the centre is ensured by 

^ > ^r- (265) 

The value of z x [xq) is positive, while b < 0. Hence, £ is given by Eq. (242), like in Ref. [ p5[ . Condition (265) is even 
more complicated than Eq. (249) and simplifications do not seem possible. This is true even for the uncharged case, 
discussed in the mentioned reference. We shall remark only that the simple case np\ = Kpo = Xo = 1, K = has 
positive pressure profile. For small enough K, the charged case must behave the same way. This ends the discussion 
of the case 1 = 3. 

Let us study next the case I = 2 when z is given by Eq. (220). Eq. (224) is not soluble unless b = —npi/5 = 
meaning that the density is constant. Then Eq. (224) becomes a particular case of the hypergeometric equation 

x (x - 1) Vxx + iivi + m + 1) x - m] y x + vmy = o, (266) 

where \ — ax , Vi + V2 + 1 = 1/2, r/3 = 0, 771772 = K 2 /2a. These relations give 



and 7712 S [—1/2, 0). The reality of rji is ensured by 

k Po > 2SK 2 . (268) 

Thus, in the charged case it is not possible to put p = po = 0. This case was studied by Wilson p3) who did not 
recognize the appearance of the hypergeometric function. He assumed that the pressure was positive and developed 
a series expansion for y. Now, since 773 = 0, the two fundamental solutions of Eq. (266) are |^9| 

Vi=xF\rh. + l,\-Vi,%x)> (269) 

2/2=FL,-i-7 ?ll i,l-x). (270) 

This is because 1 + rji + 772, ±??i,2 are not integers. The first solution is unphysical since e v vanishes at r = 0. The 
second has a finite limit 

V2(0) = Vn f rn4 _ r , (271) 

so a linear combination between the two solutions is a candidate for a regular metric. 

In Sec. VI we have shown that when e v is given by the simple ansatz (145) with k — 1, z is expressed through 
hypergeometric functions. Here we see that, in a rather symmetrical way, when z is given by the similar ansatz (50) 
(or (220) with p\ = 0), y becomes a hypergeometric function. In this paper the emphasis is laid on solutions in 
elementary functions, so we shall not pursue this issue further, except in the case where K = 0. Then Eq. (266) 
becomes 

(x-i)y xx + \y x = o (272) 
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and is easily solved 

r ... i /a 1 2 

(273) 



2 (1 - ar 2 ) 1 ' 2 + ai 



This is the expected Schwarzschild interior solution j5g] , contrary to the claims in Ref. |2J 
A model for a superdense star with the ansatz 



- = T 1 -^ ( 274 ) 

has been discussed in the uncharged case for a\ — 2 |6(]|3l|], oi = 7 |62) and a set of discrete values for a\ [ ^3| . 
Recently, it was studied for arbitrary ai both in the uncharged Q] and the charged case ]4q| . The charge function is 
chosen to be 

W (275) 



» 1 i 2 ' 

which simplifies considerably the equation for Introducing the variable 



,= [-^]'\l-a 2 r 2 ) 1/2 , (276) 



one obtains 



(! - V 2 ) Unvv - Wnn + d iVv = °> ( 277 ) 
2X 2 

di=ai+ 2 , (278) 

ai 

where d\ > is required. 

Formally, Eq. (277) is a subcase of Eq. (224) with a = 0, b = 1, x — ► r, y - * 2/77) so we can use the machinery 
developed there. Eq. (227) gives £ = arcsin?7. We may change the variable to <5 = arccos?7 because £ + <5 = 7r/2 and 
Eq. (226) is invariant under this change. Hence, y v is given by Eq. (236) and y may be found by integration. Use of 
trigonometric equalities helps to find an expression with two terms 

= / cos[( % /^+l) ( 5 + a3 ] cos [(ygl - 1) S + 03] | 

where the notation Ci = — 2a4sina3, C2 = 2a4COsa3 has been used. This formula was found in Refs. p4| , ff"5"| by 
resorting first to Gegenbauer functions and then to Chebyshev polynomials. The peculiarities of Eq. (224) permit a 
straightforward derivation in elementary functions. 

Finally, let us discuss a solution of Nduka |3f|, belonging to the class (A, q). He chose a constant z 1 and 
generalized the uncharged solution of Ref. |34j by taking 1 = 1, i.e. q = Kr. Then Eq. (24) turns into the Euler 
equation. Repeating the analysis referring to Ref. |37j, after Eq. (157), we see that y is singular at r = 0. The same 
is true for the density, p ~ 1/r 2 . 

X. THE CASE (A, TV) 

In this case the fluid satisfies the linear equation of state (17) and some ansatz for A is given. The other metric 
component is found from Eq. (27) - a linear second-order equation for y. 

Let us study first the case z = c < 1 which leads to singular pressure and density, but provides the opportunity to 
generalize the well-known KT solution |ll 42 4^j6^,Q. Eq. (27) becomes the Euler equation when po — (7-law) 



r 2 y" + k ry' - hy = 0, (280) 



3 — n , 1 — c . „ . 

k Q = ——, fci = , (281) 

71+1 C 
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whose solution is 



e v =dr 



l-fe +2fe 2 



(282) 



(1 - kof + 4fo 



1/2 



(283) 



It can be checked that 1 — ko + 2k 2 > 0. We have included only the non-divergent at r = solution. The density and 
the charge function are given by Eqs. (166), (168) respectively 



(n + 1) up = 



(l-k Q + 2k 2 )c 



(284) 



4n+(n + l) 
1 ^ — c 



(n + 1)" 

The density is positive and has a pole at the centre. Obviously 

(n + 1) 2 



c < 



4n + (n + 1) 



(285) 



(286) 



should hold. This solution is a charged generalization of the KT solution which is recovered when the equality holds. 
Due to po — it does not have a boundary and is not asymptotically flat. Another generalization was offered in Ref. 
pO] but it alters many of its properties, including the equation of state. If po + 1 Eq. (27) becomes 



r 2 y" + k ry' + (-fci + k 3 r 2 ) y = 0, 



(287) 



2np 
(n + l)c 



Its solution is given by Bessel functions 



V = r 2 



CiJ k2 (Vhr) +C 2 Y k2 (Vfc^r) 



(288) 



(289) 



Let us next discuss the ansatz (50), imposed by the condition of constant T§, like in Sec. III. Introducing x = ar 2 
turns Eq. (27) into 



x(x-l) y xx + 



n + 5 



2 (n + 1) n+1 



y x + Dy = 0, 



D= , 1 -(3n+l-^). 
2 (n + 1) V a J 



(290) 



(291) 



:(n + l) 

This is once again the hypergeometric equation and one of its fundamental solutions is y\ — F (r/i, 772 j "3; 2 1 ) where 



?72 =e-??i, 



(292) 



n + 1 



(293) 



771772 = D, 

3 — n 
= 2 (n+1)' 



(294) 
(295) 
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The solution is determined by the parameter n and the ratio Kpo/a which enters the definition of D. Eq. (294) 
shows that we can take instead rji as a second independent parameter. Eq. (291) provides one upper limit for D: 
D < aX^+IJ- In the Physical range < n < 1 it increases from 1/2 to 1. Eq. (294) is a quadratic equation for 771 and 

it is real when D < (3 — n) 2 /16 (n + l) 2 . This second limit decreases when n increases, from 9/16 to 1/16. There is 
a value no = 0.026 where both limits meet. Thus, the first inequality is stronger when n £ [0,no], while the second 
dominates for n G [no, 1]. The interval [0,no] is also the range of n where m is real when po = 0. The case n = 1 is 
degenerate because 773 = 1. This affects the second fundamental solution |59) (Vol.1). When n < 1 it is 

—1 / n - 1 1 2n \ , 

y 2=x — F [ m + —-,-- m ,——;x\ 296 
\ n + 1 2 n + 1 J 

and has a pole at r = 0. Therefore, we shall discard it. When n = 1 and rji, \ + ?7i are not integers, we have 

2/2 = FL,i- m! i;l-^ . (297) 
This function diverges at x — due to the formula 

^^^W^fraH" (298) 

r(«3-«i)r(7/ 3 -?72) 

When rji = —k or — fc + ^, fc being a positive integer, 7/2 is given by the so-called Kummer series and is reduced 
essentially to a polynomial. When 771 — k or k + | the same is true, due to the symmetry of F with respect to its 
first two arguments. In the following we shall use mainly y\ which is well-defined by the hypergeometric series and 
converges for x € [0, 1]. The generic case may be divided into three subcases: 1) D = 0, 2) D < and 3) D > 0. We 
shall search them in turn for elementary solutions. 

1) D = 0. Then rji = 0, 772 = £, npo/a — 3n + 1 and y is expressed through the incomplete /3-function 

y = C 1 + C 2 Jx-^{l-xy 1/2 dx. (299) 
When n = 1 the integral is elementary and reads 

In 1 (1 Xj 1/2 . (300) 
l + (l-a;) 1/2 

The logarithmic singularity at the beginning is obvious. When n < 1 it becomes a pole so that C2 = is required. 
Then, however, q = and the solution is uncharged. It is easy to find that up — ~a and p + ip = 0. We obtain ESU 
p0| , |52"| and Eq. (299) represents its ill-defined generalization to the charged case. 

2) D < 0. Ambiguity is fixed by accepting that 771 < 0, 7/2 > 0. A rather special subcase is 7/1 = —1/2. Then 
rj 2 = 773 and 

e" = l-a: = z. (301) 

Eq. (10) gives p + p = while 5 = 0, p = ^pj- Thus for any n we obtain the de Sitter metric whose charged 
generalizations were discussed in Sec. IV. This effective change in the equation of state was necessary because Eq. 
(27) is undefined for n — —1, and is possible when p and p are constant. It was noticed in the uncharged case for 
ESU in Ref. |47). 

Another interesting elementary subcase is given by 771 = — k when F becomes a polynomial obtained from the first 
k + 1 terms of the hypergeometric series. It may be written also as 



F(-k,e + k, m ;x) = *" ^-W/a) (1 _ 2x) = j^f^C% (VT^) , C$02. 



where (a) k = a (a + 1) ... (a + k — 1) , P and C are respectively Jacobi and Gegenbauer orthogonal polynomials |5^] 
(Vol.2). Unfortunately, this nice solution is plagued, like the de Sitter solution, and like all solutions with negative 
D, by tension. To show this let us plug our ansatz in Eqs. (166)-(168) and replace npo/a by D through Eq. (291): 

(n+l)q =2 {n+l)D-2{l-x)v x , (303) 
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(" +1 ) Kp = -( n +l) + 2(n+l)D + 2n{l-x) v x . (304) 

Since D < 0, q 2 > only if u x < 0. But then all three terms in the pressure are negative. 

3) D > 0. Then only the combination 771,2 > is possible, for when 771 < 0, 772 > 0. Eqs. (292), (294) yield 
< rji < e. The right limit decreases with n and always e < 3/2. An interesting subcase is n = 1/3 . We have e = 1, 
773 = 1/2, D = 771 (1 — 771). Then y\ degenerates into 

, = <" (1 - 2 l'- ) * . (305, 
(1 - 771) sm2i? 

x = sin 2 R. (306) 

The ranges are < 771 < 1, < D < 1/4, 4/3 < np /a < 2. When 771 becomes bigger than 1/2 the sign of C\ should 
be changed to keep y positive. Eqs. (303), (304) become 

\d-\{1-x)v x , (307) 



3ar 4 3 2 

^ = \{\-x)v x -l + 2D. (308) 

CI £ 

We demand the positivity of p. This means at least 

\{l-x)u x >\, (309) 

i.e. u x must be positive. Then, however, the first term in Eq. (307) cannot compensate the second because 2D/3 < 1/6. 
If we make q 2 positive, the pressure becomes negative in complete analogy with the case D < 0. 

Finally, let us point out that Eq. (27) may be obtained by linearization of Eq. (33). Let us search for simple 
solutions of this equation. In Sec. IV we solved the case n = — 1, when Y = using Eq. (30). Let us assume now 
that 

l(^M')' = 7 i^ r . (310) 
r 4 K ' (n+l)r 

Quite interestingly, this again leads to the ansatz (50) 

z=1 - — — — — — — -t 2 (311) 

3(77+1) V ' 

and K Po /a = 3 (n + 1), D = Then Eqs. (292), (294) give 

so that 771 and 77 are not independent. The first root leads to 771 = 773 and Eq. (301). This is the trivial Y — solution 
of Eq. (33). The other root gives 771 < and D < 0, leading to negative pressure. Something more, from the general 
solution of Eq. (33) 

Y = (313) 

z i/2 r — (C+y z-3/ 2 r^TTdr) 

it is seen that Y (and the pressure) has a pole at the centre, even when C = 0. In conclusion, we have been unable 
to find realistic solutions with positive pressure for the ansatz (50). 
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XI. DISCUSSION AND CONCLUSIONS 



Charged static perfect fluids have attracted considerably less attention than the uncharged ones, the number of 
papers being roughly an order of magnitude smaller. Solutions have been rare and the introduction hints how diverse 
were the approaches to their study. Charged dust occupies the first place in popularity and puts forth the idea 
of point-like particles. The rest is a mixture of generalizations: of the Schwarzschild idea about incompressibility, 
which may mean constant Tq or constant density in the charged case, with the limiting case of vanishing density 
and electromagnetic mass models, generalization of the idea about vacuum polarization which brings forth the de 
Sitter solution, generalization simply of well-known uncharged solutions like those of Adler, Kuchowicz, Klein, Mehra, 
Vaidya-Tikekar and others, generalization of Weyl type connections, well studied in the electrovac case. 

In this paper we have tried to show that the charged case has a life of its own when subjected to a natural and 
organized classification scheme. Surprisingly, in many respects it looks simpler than the uncharged case. The presence 
of the charge function serves as a safety valve, which absorbs much of the fine tuning, necessary in the uncharged 
case. The general formulae derived here show that the abundance of elementary solutions is probably bigger than 
in the traditional case. The proposed scheme becomes rather trivial there, representing ansatze mainly for A, v and 
sometimes p. In the charged case, however, it allows to sort out the different ideas mentioned above. Thus constant Tq 
leads to the cases (A = 1 — ar 2 , *), while constant density - to the much more difficult cases (v, p) and (A, p), similar 
to (y, n) and (A, n). Electromagnetic mass models are subcases, often spoiled by negative pressure. The point-like 
idea seems to be viable only for CD, leading to flat spacetime when the pressure does not vanish. The generalization 
of uncharged solutions firmly occupies the cases (A, q) and (z/, q). Models with n = — 1 are intimately related to the 
soluble case and are much richer than their traditional protagonist - the de Sitter solution. The other models 

with negative pressure, which seem worth being studied, are the generalizations of ESU with n = —1/3. Finally, 
the Weyl type relations, so successful in electrovac and CD environments, seem rather out of place in a 'pressurized' 
perfect fluid. 

Another advantage of this classification is that it delineates the degree of difficulty and the most convenient points 
for attack of the problem. The easiest cases seem to be (A, Y) and (u, Y). If we insist on authentic fluid characteristics 
and not their combinations, then (v,p) and (v,q) are the best candidates. On the other hand, the most difficult 
are (p, p) and (p, q) , which seem to be accessible only numerically. The most unpredictable case is probably (A, v) 
since there is no control on exactly those characteristics which must satisfy the majority of regularity and positivity 
conditions. 

A different division is between general and special cases. As was noticed in the introduction, a known solution 
belongs to any of the general cases; a solution of (v,q) may be written as (X,p) or even as (p,p) solution. Only the 
simplicity of two of the five functions v 1 A, p,p, q betrays where was the starting point. The special cases we discussed, 
apart from the ansatz (50), are essentially three: p + p = and (A, n), [v, n). All of them have a linear equation of 
state. It should be mentioned that all integrable uncharged solutions with 7- law found in Ref. [[l6| have their charged 
generalizations with p$ ^ and, hence, have a boundary. The case n = 0, po — represents CD and arises as an 
electrification of the trivial uncharged dust solution (flat spacetime) as seen from Eqs. (95)- (98). Its boundary may 
be put anywhere, since p — everywhere. The case n = — 1 generalizes the de Sitter solution into a bunch of new 
solutions, characterized by the mass function, subjected to several mild restrictions. It was completely solved in Sec. 
IV. A second, this time ill-defined generalization, is given by Eq. (313). The de Sitter solution appears also as a 
special case in Sec. X (the class D < with 771 = —1/2). Since Eq. (33) excludes the case n = — 1 this becomes 
possible thanks to an effective change in the equation of state. The option n = — 1/3 in the uncharged case is a mark 
of ESU. One possible generalization of it was performed in Sec. VII (see Eqs. (176), (177) and their discussion). The 
ESU appears also as the special class D = in Sec. X. The case n — —1/5 is special in the charged theory too, but 
its metric is singular, as was found in Eqs. (202), (203). The KT solution is generalized in Eqs. (282)-(286). In the 
uncharged case these exhaust the integrable cases of a complicated Abel differential equation of the second kind. In 
the charged case it is replaced by a linear first-order Eq. (19), which has many other solutions, discussed in Sec. VII, 
or by a linear second-order Eq. (27), discussed in Sec. X. 

From a mathematical viewpoint, the charged case delivers a surprising variety of equations like those of Euler (157), 
(280), Bernoulli (75), (313), Riccati (30)-(33), Emden-Fowler (113), Abel (37), (102), (114) and the hypergeometric 
equation (266), (290). The cases containing A quickly lead to special functions whose spectrum is also rich: Ei, en, 
Bessel, hypergeometric, incomplete /3-functions, Jacobi, Gegenbauer and Chebyshev polynomials, etc. Therefore, a 
more thorough study would require computer simulations of special functions, which are easier than purely numeric 
simulations. 

From a physical viewpoint, the most interesting results are, in our opinion, the following: 

1) The case of constant Tq is solved by a simple algorithm involving algebraic operations, one differentiation and 
one simple integration upon a generating function Y , which satisfies few simple inequalities, Eqs. (51)-(53). 
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2) The charged de Sitter case is completely soluble in terms of the mass or the charge function, Eqs. (56)- (59). An 
example is given by Eqs. (71)-(74). There is a generalization with positive pressure in Eqs. (76)-(78). An example is 
given by Eqs. (79)-(82). 

3) The general solution for a given pressure may be written as three contributions to z; from regular CD, from 
general CD and from the pressure, see Eq. (83). There is a general CD solution which creates a halo around charged 
fluid balls and may postpone their junction to a RN solution. An example is given by Eqs. (115)-(140). 

4) The solutions of Korkina-Durgapal possess realistic charge generalizations given by Eqs. (142)-(145). An example 
is given by Eqs. (148)-(153). 

5) Fluids with linear equation of state are integrable for any n as seen from Eqs. (19), (27), (33), (159)-(165). 
Elementary solutions may be found by several simplification techniques. Physically realistic is an infinite series of 
models given by Eqs. (173), (174), (178)-(180). Realistic examples are given by Eqs. (182), (185). The status of the 
model given by Eq. (190) is still unknown. 

6) The case (y,p) is closely connected to {v,n) and is much more difficult than its companion with constant Tq, 
even when the density is constant or zero. An example with special functions is given by Eq. (213), another one - by 
Eq. (216). Their physical status is still undetermined. 

7) The (p, q) case discussed in Ref. |^6| is incomplete and has errors. It is treated in detail in Sec. IX and its 
connection with Ref. 0] is elucidated. It has several realistic cases, but electromagnetic mass models all seem to have 
negative pressure. Realistic examples are given by Eqs. (218), (234), (238)-(245) and Eqs. (263)-(265) with the values 
specified after them. A formal, but intriguing parallel is drawn to star models with simple spatial geometry |45|. 

8) The solution of Wilson Q is expressed by a hypergeometric function and reduces to the Schwarzschild interior 
solution when the charge vanishes. Its realistic status is not known. 

9) The KT solution has a charged generalization which is likewise singular, but may serve as a focal point to some 
regular solution, see Eqs. (282)- (286). 

10) Solutions of the case (A, n) with constant Tq are expressed in hypergeometric functions, see Eqs. (290), (291). 
All degenerate elementary solutions have either negative pressure or singularities. 
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